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Abstract 

The fully nonlinear development of small wavelength Gortler vortices in a growing 
boundary layer is investigated using a combination of asymptotic and numerical methods. 
The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who 
discussed the initial development of small amplitude vortices in a neighbourhood of the 
location where they first become linearly unstable. That development is unusual in the 
context of nonlinear stability theory in that it is not described by the Stuart-Watson 
approach. In fact the development is governed by a pair of coupled nonlinear partial 
differential evolution equations for the vortex flow and the mean flow correction. Here the 
further development of this interaction is considered for vortices so large that the mean flow 
correction driven by them is as large as the basic state. Surprisingly it is found that such 
a nonlinear interaction can still be described by asymptotic means. It is shown that the 
vortices spread out across the boundary layer and effectively drive the boundary layer. In 
fact the system obtained by writing down the equations for the fundamental component of 
the vortex generate a differential equation for the basic state. Thus the mean flow adjusts 
so as to make these large amplitude vortices locally neutral. Moreover in the region where 
the vortices exist the mean flow has a ‘square-root’ profile and the vortex velocity field 
can be written down in closed form. The upper and lower boundaries of the region of 
vortex activity are determined by a free-boundary problem involving the boundary layer 
equations. In general it is found that this region untimately includes almost all of the 
original boundary layer and much of the free-stream. In this situation the mean flow has 
essentially no relationship to the flow which exists in the absence of the vortices. 


This work was supported under the National Aeronautics and Space Administration 
under NASA Contract No. NAS1-18107 while the authors were in residence at ICASE , 
NASA Langley Research Center , Hampton ,VA 23665. 
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Introduction 

Our concern is with large amplitude Gortler vortices in viscous incompressible flows 
over walls of variable curvature. Much of the recent interest in the Gortler instability 
has been motivated by practical problems such as Laminar Flow Control or the flow over 
turbine blades. A particular cause for concern in Laminar Flow Control is the question 
of whether Gortler vortices are likely to induce premature transition because of their 
effect on the receptivity of the original boundary layer to Tollmien-Schlichting waves or 
crossflow vortices. As yet little progress has been made with the latter problem but 
the corresponding internal fully developed flow problem has been discussed by Hall and 
Bennett (1986), Bennett and Hall, (1987) and Hall and Smith (1987). Thus, for example, 
Bennett and Hall (1987) showed that fully, nonlinear Taylor vortices in a curved channel 
can have a massive influence on the growth of lower branch Tollmien-Schlichting waves. 
The absence of a fully nonlinear description of Gortler vortices in external flows means that 
a similar interaction problem cannot yet be attempted for these flows. Here we shall show 
how the required large amplitude vortex flows can be described asymptotically. However 
it is perhaps useful if we first describe earlier work on the Gortler problem. 

The original work of Gortler (1940) and the later work of, for example, Hammerlin 
(1955, 1956), Smith (1955) and Floryan and Saric (1979) did not take care of boundary 
layer growth in a selfconsistent manner. Unlike the corresponding Tollmien-Schlichting 
instability problem, where the fact that instability occurs at relatively high Reynolds 
numbers renders parallel flow calculations reasonably accurate, the neglect of non-parallel 
effects in the Gortler problem leads to inconsistent and, in some cases, physically absurd 
results. Thus as extreme examples some of the results predicted by these calculations 
showed instability at zero Gortler number or zero wavenumber. 

In fact at high wavenumbers the above theories gave consistent results and Hall (1982a) 
showed that in this regime an asymptotic solution of the non-parallel problem is possible. 
The vortices are found to become linearly unstable at a particular downstream location and 
to concentrate themselves at some depth in the boundary layer. This depth corresponds to 
the position where Rayleigh’s inviscid instability criterion for the boundary layer is most 
violated. At O(l) wavenumbers no asymptotic or self-consistent parallel flow calculation is 
possible and the linear instability partial differential equations must be solved numerically 
as was done by Hall (1983). These equations are parabolic in the downstream variable 
and so can be solved by a marching procedure and the neutral location for a disturbance 
inserted at a given position can be calculated. Hall (1983) showed that the neutral curve 
produced by such a calculation depends on the nature of the initial disturbance. However 
each such neutral curve has none of the obvious anomalies of the parallel flow neutral 
curves and at high wavenumbers the different neutral curves merge into the asymptotic 
result. 

The initial nonlinear development of Gortler vortices of small wavelength was discussed 
by Hall (1982b). It was found that the nonlinear interaction is dominated by a ‘mean-field’ 
type theory in which the fundamental mode and mean flow correction reinforce each other. 
The higher harmonics play no role in this interaction so the Stuart-Watson approach is 
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not applicable. In order to motivate the asymptotics to be used for much larger vortices 
it is instructive to summarize the essential details of the calculation of Hall (1982b) 

Suppose then that a vortex of nondimensional wavenumber e -1 is locally neutrally 
stable at the downstream location x = x of a two-dimensional boundary layer flow (u, v). 
If the most unstable part of the boundary layer is at a depth y = y it is necessary to define 
variables X and £ by 


e 



(1.1a, b) 


and the total downstream velocity component expands in the neighbourhood of (x, y) as 


u = u(x,y)+eXu x (x,y) + e^u y (x,y) H b 

+ e*[U 1 (X,<;)exp(j) + C.C.} + --- . ^ 

For a boundary layer which has a local Gortler number increasing faster than the fourth 
power of the local wavenumber the evolution equation for u M , the mean flow correction, 
and the fundamental U\ take the form 


( * 

d > 

d 1 

w 

dX) 

UM = a? 1 


d 2 
d ^ 


2 d 1 2 „ 

3dX + 


}*" 


2 Ui 


du 


M 


dS 


(1.3a, 6) 


The linearized form of (1.3) shows that U\ can be expressed in terms of parabolic cylinder 
functions. Otherwise (1.3) must be solved numerically by a marching procedure. However 
at large values of X it was shown by Hall (1982b) that um and U\ develop a surprisingly 
simple asymptotic structure. In fact um and Ui can be written down explicitly in terms of 
X and a similarly variable £ in a finite interval £i < £ < £ 2 . Near £ 1 , £2 the fundamental 
satisfies a nonlinear Airy equation and is reduced to zero exponentially whilst um persists 
above £2 and below £ 1 . Of crucial importance is the fact that when this structure develops 
um ~ X 3 / 2 , U\ ~ X 1 / 2 so that the total mean flow correction c’^um is comparable with 
the basic state for X ~ l/e. Thus if x — x is 0(1) the vortices will drive a mean flow 
correction as big as the original mean state. 

The remarkable feature of the latter regime is that the nonlinear Gortler equations can 
still be solved asymptotically in this region. Thus we are able to describe asymptotically 
a large amplitude disturbance capable of altering at zeroth order the basic state. The 
structure of the large X solution of Hall (1982b) continues to describe the flow but with 
the major change that the depth of the fluid where vortex activity persists is now 0(1). In 
this layer the mean flow is driven by the vortices and indeed is determined as a solvability 
condition on the equations for the fundamental there. The downstream mean velocity 
component in this layer then has a simple square root profile and the mean equations 
drive a finite amplitude vortex. Thus there is an exact reversal of the usual roles of 
the equations for the mean and fundamental obtained by Fourier- analysing the spanwise 
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dependence of the disturbance. The situation is not unlike the scenario postulated some 
years ago for turbulent flows by Malkus (1956). We recall that Malkus argued that the 
‘mean’ part of a turbulent flow would organize itself so that any ‘modes’ were marginally 
stable. We shall refer to the layer where the vortex activity is concentrated as I. This layer 
is bounded by ‘transition’ layers of depth e 2 / 3 where the vortices are again reduced to zero 
as solutions of a nonlinear Airy equation. We refer to these layers as IIa,b respectively 
and the remainder of the flow is denoted by regions IIIa,b as shown in Figure 1. 

In regions IIIa,b there is no spanwise dependence for the flow and the velocity field 
satisfies the boundary layer equations. The solutions of these equations must be matched 
with those emanating from Ila, b whose positions are unknown functions of the downstream 
variables. Thus a complete description of the flow requires the numerical solution of a free- 
boundary partial differential system. 

The region I of vortex activity in general grows as the flow moves downstream until 
it eventually occupies almost all of the extent of the boundary layer which would exist in 
the absence of the vortices. In fact the only part free of vortices is a thin layer at the 
wall and the vortices even take over part of the free stream. At this stage the flow has 
effectively no relationship with the flow which would exist in the absence of the vortices. 
The procedure adopted in the rest of this paper is as follows: in §2 the nonlinear Gortler 
vortex equations are derived. In §3 an asymptotic solution of these equations for large 
amplitude small wavelength vortices is given. In $4 an asymptotic solution of the initial 
and ultimate downstream stages of the free-boundary problem obtained in §3 is discussed. 
In §5 a numerical scheme which we have used to solve this free boundary problem is 
described. Finally in §6 we discuss our results and draw some conclusions. 

2. Formulation of the nonlinear Gortler equations. 

We consider the flow of a viscous fluid of kinematic viscosity v, density p , over a wall 
of variable concave curvature a -1 x(X/L). Here X denotes distance along the wall, a is a 
typical radius of curvature and L is a typical length scale along the wall. If Uo is a typical 
flow velocity we define a Reynolds number Re by 


Re = (2.1) 

and a curvature parameter 6 by 

6 = -. ( 2 . 2 ) 

a 

We confine our attention to the limit Re —> <x> with the Gortler number G defined by 

G = 2Rl6 (2.3) 

held fixed. We take (X, Y, Z ) to be co-ordinates along the wall, normal to the wall and in 
the spanwise direction respectively. If (U,V,W) denotes the corresponding velocity vector 
we define dimensionless co-ordinates ( x,y,z ) and velocity ( u,v,w ) by 
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and 


(U,V,W,) = Uo(u,vRe 1/ \wRe 1/: ‘). 


(2.4) 


We restrict our analysis to flows with u — ► 1, y — * oo so that the pressure P can be written 
in the form 

Ke 

The Navier-Stokes equations for the flow can be written in the form 

du dv dw 

dx dy dz 

du du du 

Vu = u— + v— + w — , 
dx dy dz 

_ 1 ~ 2 dp dv dv dv 

v ”“ 2 Gxu -T y =u ai + v ai JrW Yz' 

_ dp dw dw dw 

vw — — — = u— 1- u— b W——, 

dz dx dy dz 


(2.5a,b,c,d) 


q 2 q2 

where V = + — - and terms of relative order R Z 1 ^ 2 have been neglected. We further 

dy z dz 1 

note that we have assumed that the flow is steady. 

In the absence of any Gortler vortex we can write (u, v,w) = (tt,u,0) in which case 
u,v satisfy 

du du d 2 u 

U dx V dy dy 2 ’ 


du dv 
dx dy O’ 

u = v — 0, y = 0, u — > 1, 


( 2 . 6 ) 


v —* oo. 


u = 1, x = 0, y > 0. 


Thus (u v 0) is just a Blasius boundary layer until a Gortler vortex begins to grow at some 
downstream location 

Finally in this section we note that our analysis can easily be modified to take care of 
pressure gradient driven flows. In that case P in (2.4) must be altered so as to contain an x 
dependent component of size pUfi and the boundary conditions in (2.6) must be modified. 

We shall obtain a solution of (2.5) which satisfies 


u = v = w = 0, y = 0. 
u — > 1, y — ► oo. 

In fact we shall see that when y —*■ oo,v —*■ v[x),w — ► 0; in general v(x) is not the 
y velocity component at infinity of a Blasius boundary layer so the higher order outer 
‘inviscid’ problem associated with (2.5) is modified by the instability. 


(2.7a, b) 
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3. The evolution equations for large amplitude Gortler vortices 

We now develop an asymptotic solution of (2.5) valid in the limit of small vortex 
wavelength. We suppose that the boundary layer (u(x,y),v(x, y)) becomes linearly un- 
stable to Gortler vortices at x = x*. Furthermore we assume that the curvature function 
x(x) is such that the flow becomes more unstable with increasing x. Thus if the boundary 
layer is a Blasius boundary layer we require that x increases more rapidly than 12 . This 
restriction is discussed by Hall (1983) and is a direct consequence of the scaling of the right 
hand branch of the neutral curve for Gortler vortices. 

The discussion of §2 suggests that beyond x — x* the flow will support fully nonlinear 
Gortler vortices. The flowfield is therefofe split up as shown in Figure 1. The vortex 
activity is therefore confined between yi{x) and y^x). The layers denoted by regions 
IIa,b are required in order to smooth out the algebraically decaying vortices in region I. 
Later we will see that these layers must be of thickness e 2 / 3 in the limit e — > 0. We expand 
the Gortler number G in the form 


G = Gq6~ 4 + Gie~ 3 + • • • , 


(3.1) 


whilst in region I the appropriate expansions of u,v,w and p are 

u = u 0 + eu! + • • • + {eE{U£ + tU\ + •••) + e 2 E 2 (U 2 + eU 2 + •••)••• + C.C.} 

v — vo + ev\ + • • • + + • • •) + £®E 2 (V q + eV^ 2 + •••)••• + C.C.} 

w = E(Wq + eWl + •••)+ eE 2 (W$ + eW? + •••)••• + C.C. 

P = Po + epi + • • • + {e 1 E(Pq + ePi + • • •) + e°E 2 (V 2 + eV 2 + ••■)••• + C.C.} 

(3.2a, 6, c,d) 

Here ‘C.C.’ denotes ‘complex conjugate’ whilst E = exp (tz/e). In the absence of a vortex 
the flow reduces to the basic state w = p = 0, u = u,v = v. We note that the only z 
dependence in (3.2) is through E so that uo,vo, etc. are functions of x and y only. 

The expansions (3.1), (3.2) are then substituted into (2.5) and like powers of e for 
each Fourier component are equated. The mean flow (tZo^o) is determined from the 
zeroth order continuity and x momentum equations. The pressure field po is then found 
from the normal momentum equation. The equations to determine Uo and vo are 

du 0 dvo_ _ 0 
dx + dy ~ ’ 


duo _ duo 
Uo -s 1- Vo 


dx 


dy dy 2 


d 2 u 0 = v ,d_ui + %l eui_ Woiij x + 


dy 


dy 


(3.3a, b) 
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Thus the boundary layer equations are now driven by the nonlinear interaction of the 
Gortler vortex with itself. The zeroth order equations for the fundamental terms in the 
core are 

dV' 


dy 


0- + iW^ = 0, 


U i +F i^-o 
0 + 0 dy ’ 

Vo 1 + GoX^o^o = 

-iPo = Wo- 

The consistency of equations (3.4b,c) requires that throughout the core 

G oX«o-5— = 

dy 


(3.4a, b, c, d ) 


(3.5) 


and (3.3) then becomes 


dup dvp 
dx dy 


Up 


du o 
dx 


+ V 0 


du 0 

dy 


d 2 u 0 

~dy r 


a auo 

dy X dy 


Vp 1 12 


}• 


(3.6a, 6) 


Thus (3.5) determines the mean flow streamwise velocity component in the core whilst 
(3.6b) then determines | Vq 1 | 2 . It follows that in the core the boundary layer flow is now 
being forced by the vortex which from (3.6b) is itself driven by the boundary layer. This 
is, of course, the exact reverse of the roles played by these equations in the linear or weakly 
nonlinear descriptions of this problem. In fact the interaction described by (3.5) and (3.6) 
is not unlike that postulated some years ago by Malkus (1956) for turbulent flows. In the 
latter theory it was argued that the ’mean flow’ for a turbulent flow organizes itself so as 
to make any possible instability marginally stable. In the large amplitude Gortler problem 
this mean flow adjustment is achieved by centrifugal effects and is described by (3.5). We 
can integrate (3.5), (3.6) to give 


y/ a(z) + 2y 
U °“ V^oX 


(3.7a, b) 


-a'(x)s/a(x) +2y {a(z)+2y}lx' u _ ^ 

Vo — + <tyUoxx — {) 


where a(x) and b(x) are arbitrary functions of x and a dash denotes a derivative with 
respect to x. The function | Vq | 2 is found by integrating (3.6) to give 


B{x) + 2 


dup 

dy 




by/d + 2 y + 1 (a + 2y) 2 x / 4 1 

/GpX 12 Gpx 2 y/a + 2 y/ GoX 


(3.8) 
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where B{x ) is another function of x to be determined. The function | Vq | 2 cannot be 
negative so 3/1 and y 2 which determine the edges of I satisfy (3.8) with Vq = 0. If B(x) is 
then eliminated we obtain 

by/a + 2yi + 1 (a + 2yi)V + 1 = 

V G oX 12 G 0 X 2 v /a + 2yiV /G ’oX ^ 

fry a + 2y 2 | 1 (q + 2y 2 ) 2 x' | 1 

\/GoX 12 C0X 2 y« + 2y 2 \/G ! oX 

The above equation is not of course sufficient to determine a, b, 3/1 and y 2 so we are not yet 
able to determine the location of the layers Ila, b. The thickness of these layers is deter- 
mined by a balance between diffusion across the layers and convection in the streamwise 
direction. This balance shows that the layers are of thickness e 2 / 3 so that in Ila we write 


5 ” e 2/3 ' 


d d d 

Hence in Ila we replace — and — by — — 

ox ay ox 

12 


y 2 


and 


1 d 


respectively. We can 


e 2 / 3 di e 2/3 ae 

see from (3.8) that | V/ | 2 ~ 3/2 — 3/ when y — * 3/2- in which case the fundamental velocity 
components and pressure in (3.2) decrease by 0(e 1 / 3 ) where 3/ — 3/2 = 0(e 2 / 3 ). The 
appropriate expansions in Ila therefore take the form 


u — uo + e 2 / 3 ui + • • ’-\-{t 4 ^ 3 E{JJq\ -\-e 2 ^ 3 Un + • • •) -|- e 8 ^ 3 JE 7 2 (U02 + e 2 ^ 3 C 12 + * ••)+•• • + C.C.} 

v = Vo + e 2 ^ 3 u 1 + • • * + {c 2 / 3 E(V 0 i+e 2 / 3 Vu + - • •) +e 2 ^ 3 jE/ 2 (V r 02 + e 2 ^ 3 k r i 2 + ,,- ) + -,- + C.C.} 

w = e~ 1 / 3 E(Woi + c 2 / 3 W n + •••)+ e 1/3 E 2 {W 02 + e 2 / 3 W 12 + •••)+••• + C.C.} 

P = p 0 + e 2,3 p x H 1- {e~ 4 / 3 E(P 0 i + e 2 ^ 3 Pu + • • •) (3.10) 

+ e 2 / 3 E 2 {Po2 + f 2 / 3 T > i2 + •••) + •■• + C.C.} 

The coefficient in the above expansions are functions of x and £. 

It has been anticipated above that the first harmonic functions also decrease in size 
in Ila ; this decrease is forced by the form of the equations for U 2 , V 2 in the core. The 
appropriate equations are 

4U 2 + v ^ = F u (3.11a, b) 

4V 0 2 + G oX U 2 u 0 = F 2 , 

where F\ and F 2 are quadratic in Uq, Vq. The forcing terms are therefore 0(e 2 / 3 ) in Ila 
and so Uq, V 2 , Wq, Pq must be rescaled in that layer as indicated in (3.10). We note that 
(3.11) can always be solved in the core since Uq satisfies (3.5). A similar analysis applies 
to the higher harmonics so we conclude that the fundamental effectively drives the mean 
flow and all the higher harmonics. The situation is quite different from the usual type 
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of weakly nonlinear based on the Stuart-Watson method where the mean flow remains 
essentially identical to that present before any instability occurs. 

We now return to the solution of (2.3) in the transition layer Ila. The first two terms 
in the expansion of the mean flow in this layer satisfy 


d 2 u 0 

dp 


= 0, 


d 2 u i 
dp 


= 0 


and the solutions of these equations which match with the solution in the core are 


_ y/a + 2 y 2 

U °' VGoX 


u 1 


£ 

V^oxV^ + 2y 2 


(3.12a, b) 


The function u 2 is forced by the fundamental terms in (3.10) so that we cannot proceed 
further without solving for the latter terms. The equations to determine (Ucn, Voi) and 
(^lijk'u) are found to be 


Ur 


+Voi 


01 • 


= 0 , 


VGoX-s/a + 2y 2 
Vox + Uqi\/ Gqx^J a + 2t/ 2 = 0, 
Vn 


u 11 + 


V^oxV® + 2y2 


T r ^ u 2 , d 2 U 01 
01 ae + dp 


(3.13a, b, c, d) 


Vix + Uii\/G^xy/aT2y^ = U 0l + 2^i 

y/a + 2 y 2 dp 


The equations (3.13a,b) are of course always consistent but (3.13c,d) are only consistent 
if an orthogonality condition holds. To obtain this condition we must proceed to higher 
order in the mean flow expansion. 

The first order solution for vb is found to be 


-a/y/a + 2 y 2 {a + 2y 2 } 3/2 x / __ b 

2\/^oX b\/ GqXX 


whilst u 2 satisfies 


a 2 u 2 _ ~x> 2 + 2y 2 ] 6 2 d_ 

%X 2 Gq y/Go X\/a + 2 y 2 ) y/^oXy/a + 2y 2 a£ 

which can be integrated once to give 


du 2 

~di 


= -t 


x'\ a + 2y 2 ] 


+ 


3x 2 ^o yj GqX\/ a + 2y 2 


2 I Vr 


01 


V GoxV a + 2y 2 


+ /• 


(3.14) 
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Here f(x) is another function which can only be determined at higher order; for our 
purposes here it is not required here. Having determined — ■£. we can now write down the 


appropriate solvability condition for (3.13b): 




+ S(x) £V 01 = | | Voi | 2 Voi V^+2^y/G^ + G oX {a + 2y 2 ) V 01 /. (3.15) 


where 

5 (*) = a + 2v ~ + b VGoxV a + 2y 2 }. 

This equation is a particular form of the second Painleve transcendent and has been shown 
by Hastings and Mcleod (1978) to have a solution such that 

~^\Vqi\ 2 V^+^V^oX* £ -°o, | Voi HO, 


It follows that in Ila the fundamental terms decay to zero so that the finite amplitude 
Gortler vortex is trapped below region Ilia. We note that a similar analysis for the higher 
harmonics shows that these functions also decay exponentially to zero in Ila. However 
the mean flow is virtually unaltered by the presence of Ila, thus the first two terms in the 
expansion of the mean flow in Ila are simply obtained by expanding the mean flow in I 
in terms of £. This means that the mean flow in Ilia must to zeroth order have u,u y and 
v defined by the coreflow solution evaluated with y = j/2- 

An identical analysis to that above shows that the 2 -dependent part of the flow is 
reduced to zero exponentially in lib. Hence in Illb there is only a mean velocity field. 
Thus in IIIa,b we write 


u = u + 0(e 2//3 ), v = v + 0(e 2 / 3 ), w = p= 0, 

so that (u,u) defined in (0, yi) and (y 2 ,oo) satisfies 

_du _du d 2 u 

U dx V dy dy 2 ’ 

du dv 
dx "I" dy 

These equations must be solved subject to 

u = v = 0, y — 0, 


and 


u-> 1, y -* oo, 


yfa + 2 yy _ _ 1 

VGoX 5 y y/GoXy/ a + 2 yy 5 


(3.16) 


(3.17) 


(3.18) 

(3.19) 


(3.20) 
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g / v /o + 2yy | {g + 2y J } 3 / 2 x / 
2\/G 0 X 6\A?oXX 


y = yj, i = 1 , 2 . 


(3.21) 


The equations (3.17) - (3.21) together with (3.9) specify a free boundary problem for yi,y 2 
and the functions a(x ) ,6(x) . Clearly no analytical solution to this problem is available. In 
fact it is also necessary to specify something about the ‘upstream’ nature of the instability. 
This amounts to finding an asymptotic form for the solution of this system close to the 
value of x where the original boundary layer becomes linearly unstable. This will also be 
discussed in §5 where a large x solution of the system will also be developed. In §5 a scheme 
which we have used to solve the system numerically will be described. However there is a 
special case where a similarity solution of the system can be generated. This corresponds 
to the case when the undisturbed boundary layer is Blasius flow and x is proportional to 
x 1 / 2 . Though this particular case is perhaps not of much practical importance its solution 
is instructive in that it suggests how to solve the full free-boundary partial differential 
system numerically. For that reason we now indicate how the similarity solution can be 
calculated. 


Consider then the curvature distribution x{ x ) defined by 


x(x) = \pLx. 


(3.22) 


If the undisturbed flow is a Blasius boundary layer the local Gortler varies with x 2 
as does the fourth power of the local spanwise wavenumber. This is consistent with the 
scaling of the right hand branch of the neutral curve discussed by Hall (1982a). Hence on 
the basis of linear theory the flow is either stable or unstable for all values of x. We define 
the similarity variable r\ by 



(3.23) 


and seek a solution of (3.17) by writing 


« = /'(» 7 ), V= Vfx^ /_ ^' 


(3.24) 


The functions a and b then take the form 


a — aV2x, b 


b 

~\/2x’ 


(3.25) 


We can then show that the free boundary problem specified by (3.9), (3.17) and (3.21) is 
equivalent to 

/"' + = 0> (3.26) 


on 


0 < rj < rji,r )2 < rj < oo 


subject to 


f = f = 0 , V = 0 , 


(3.27) 
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and 

f'(oo) = 1, (3.28) 

together with the matching conditions 


f'M 


yfa-\-2rj ~j 

'Wo 


fivj) = Vjf'iVj) - 


a>/a + 2r7y (q + 277 J -) 3 / 2 ~ 

2y/G^0 + 6 \/Uq 


and the scaled form of (3.9): 


J = 1,2, (3.29a, b) 


by/a + 2rjl + ^ ^ {a + 2r?i} 2 + 


\fG o 


Gc 


1 


y/a + 2rfiy/Go 

6y/d + 2 T ?2 1 {a + 2 t7 2 } 2 1 3Q ^ 

y/Go 12 Go 


Here we have replaced — /= by rjj for j = 1,2. Thus in the special case x = V^x we have 

V2x 

a free boundary ordinary differential system to solve. The unknowns are a,b and 
the boundaries of the region where a finite amplitude vortex exists. These constants and 
the function / were obtained by following the procedure shown below. 


(i) Integrate (3.26) forward from rj = 0 with /( 0) = /'(O) = 0 and an initial guess for 
/"( 0). This integration is stopped at 77 = r?i where / , (f?i)/ ,, (»7i) = 


(ii) Using calculate a, b from (3.29a,b) with j = 1. 

(iii) Now calculate rj 2 from (3.30). 

(iv) Finally integrate (3.26) from 77 = r} 2 to some suitably large value of rj with 

f(r ) 2 ), defined by (3.29a, b) with j = 2 and using 


/'(*?2)/"(»?2) 



(v) If f(oo) / 1 we return to step (i) and alter the initial value of f" at the wall in 
some suitable manner. 

The above procedure was found to converge when used with Newton’s method to 
update the value of f" at the wall. We found that, as expected, a solution exists only for 
G 0 > 4.2 the neutral value of G 0 . Beyond G 0 = 4.2 the values of rj 1 and r) 2 respectively 
decrease and increase from their limiting value »h = r) 2 = 1.56 when Go = 4.2. At larger 
values of Go, rj 2 varies linearly with Go whilst rji goes to zero like Gq 3 . These scalings 
can be recovered by an asymptotic investigation of (3.26) — (3.30) in the limit Gq 00 . 
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Such a calculation shows that for Go » 1, 772 — |G 0 , rj 1 ~ .12 Gg 3 and the value of 

772 predicted in this limit is shown in Figure 2 where we have shown the result of a full 
numerical solution of (3.26) - (3.30). 

Thus for the particular case x = the partial differential system governing the 

development of the mean flow can be reduced to an ordinary differential one. Moreover the 
results shown in Figure 2 suggest the type of behaviour for yi,y 2 which we should expect 
for the general problem when x increases beyond its neutral value. At large values of x 
we shall see that a structure essentially identical to that shown in Figure 2 for Go >> 1 
is set up. Finally we note that for other similarity boundary layers we can always choose 
a particular curvature distribution x(^) which leads to an ordinary differential system for 
the mean flow. 

4. The initial and ultimate downstream developments of the instability. 

We first describe the development of the instability close to the downstream location 
where the boundary layer (u,v,0) becomes linearly unstable. The description of the flow 
in this region is similar to that given by Hall (1982b) in the weakly nonlinear regions so 
we shall not give all the details here. Suppose then that vortex with wavenumber e -1 
is linearly unstable for x > x* and the instability originates in a layer of thickness e 1 / 2 
centred on y = y*. The weakly nonlinear theory of Hall (1982b) suggests that the finite 
amplitude vortex occupies a region of depth 0(x — x*) 1 / 2 for | x — x* |<< 1. We therefore 
expand y x and y 2 in the form 

Vi — V* — {x — x*)^y + 0(x - x*), 

V2 = y* + {x - x*)?y + 0(x - x*), (4.1a, 6) 

It is convenient for us to now define a similarly variable £ by 

* _ y~y* ~ JL ( 4 21 

(X-X*)V2 “ *1/2 ^ 


so that 


d 

dx 


a 

ax 


e d 


a_ 

dy 


2 Xdi' dy X l / 2 d£‘ 
In region I we express (u 0 ,v 0 ) which satisfies (3.5) in the form 

uq = u + X z ! 2 um (£) + • • • 


v 0 = v + Xv m (£) H . (4.3a, b) 

Here (a, v, 0) is again the boundary layer flow which exists in the absence of the instability. 
In the neighbourhood of (x*,y*) u,v, and the curvature x expand as 


w — uoo + X 2 0 + X{£ 2 ii2o + «oi} + X 2 £ttn + • • • , 

v — vqo + X? £i>io + AT{£ 2 u 2 o + u oi} + + •••, 
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X — Xo + Xxi + X 2 X2 + • • • • 


(4.4a, 6, c) 


If we substitute (uo,vo) from (4.3) into (3.5) and use (4.4) then, equating like powers of 
X°, X?, X we obtain 

GoXqUqoUio = 1 , 


— Xi 

^oXo^oo 


M io "b 2 «oo «20 = 0, 


du M 

dt 


+ u n + 


3£ 2 «l0 ll 20 

u 00 


+ 


Uqo 


(4.5a, 6, c) 


Equation (4.5a) is satisfied if x* is the neutral location whilst (4.5b) holds because y* in 
the linear theory is chosen so that 


d . du . 

= °’ 


y = y 


Equation (4.5c) can be integrated once to give 

Ul Vf = A 0 f — Ai £ 3 


(4.6) 


where 


Aq — u 1 1 — 


Xi«io 

Xo 


u 01«1Q 
«00 


and 


«10«20 
Ai — 

u OO 

We have anticipated in (4.7) that u m is an odd function of f; 
been determined at a later stage. 


(4.7a, 6). 


this would otherwise have 


The Gortler vortex function 


Vq 1 in region I expands as 

I Vo 1 | 2 = »■(?)+••- 


and F(£) can then be found from the zeroth order approximation to (3.6b). This equation 
can be integrated once and F(£) is found to vanish at £ = where 


I t± | 2 = 


C + Ap 

3Ai + |-uooAo 


(4.8) 


Here C is an unknown constant of integration to be determined later. It follows from (4.1) 
and (4.8) that 

= -y, £+ = y- 

In region Ila the forcing of the mean flow ceases and [ Vq j is reduced to zero exponentially 
in the manner described in §3. It remains for us to determine the constant C appearing 
in (4.8). We can write down the mean flow in regions IIIa,b in the form 


u ~ u + X5Uas{£) + • • • 
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C = F+XVm(0+*** (4.9) 

It is easy to show that Um(£) and Vm(£) are odd and even functions of £ respectively so 
it is enough for us to consider Ilia where £ > £ 4 . . The functions U^j and Vm satisfy the 
equations 

d 2 U M 1 f . d 1 „ 


^^ + \ Uoo { ( 7(~ 3 } Um ~ 0 

0 

and in order to match with the core solution (4.6) Um must satisfy 

Um = Aof+ — Aif®, U' M — A 0 — 3Aif+, £ = £+• 


(4.9a, b) 


(4.10) 


Since Um must also tend to zero when f —* 00 (4.9a), (4.10) constitute a nonlinear eigen- 
value problem for After some manipulation we find that 


U M = 


(A 0 e + - ex P {-^[e 2 - qmb y^e) 


5 I* 

Vm = £Um — — J 


U M d<i, 


(4.11a, b) 


where U(a,x) is a parabolic cylinder function. The eigenrelation for £ 4 - is then found to 
be 


4P(|.^f+) , /u 5 > 


+ \ -?■(+> = 


Ao — 3Ai£+ 
A 0 e+-Ai4 


(4.12) 


V 2 {U&VWZ+) V 2 Ao^-A^r 

and C is then given by (4.8). Moreover t/i and 3/2 are then known correct to order (x— a;*) 1 / 2 

The functions a(x) and b(x) are found by comparing (3.7) and (4.3) with u,v,x given 
by (4.4). Thus, for example, we find that a(x) is given by 


a(x) = G 0 x |«oo - + 2Xu 00 [uoi - A 0 y * -u u ] H | (4.13) 

for X 1 . Thus we can construct an asymptotic solution of the free boundary partial 
differential system of §3. This asymptotic solution can of course be used to start off a 
numerical solution of this system by a marching procedure. 

Next we suppose that the curvature distribution function becomes large when x — * 00 . 
More precisely we suppose x( a: ) / x? — * 00 , x — >■ 00 so that the flow is unstable on the basis of 
linear theory when x — * 00 . In this limit it is convenient to perform an asymptotic analysis 
directly on the solutions obtained in §3. For definiteness we assume that x = x M ,M > | 
when x — ► 00 . 


15 


The scaling for x ~ * 00 is suggested by the large Go solution discussed in §3 and can 
be inferred from (3.7). It is easy to show that a solution of the ‘upper problem’ in Ilia 
requires that 3/2 >> x * and that Uo = O(l) when y = 3 / 2 - Thus I must be of depth O(x); 
the core solution for y — O(x) can then be written 


I 2y 

“° “ V + 


( 4 . 14 ) 


and if we take 


G J 
y 2 = -j- + 


uo — ► 1 H when y — ► 3 / 2 —. We then find that 

G 0 Mx m ~ 1 

v 0 -*• + • • • , 

when y — > 3 / 2 — .Thus a solution of (3.17) for (u,v) in the upper layer Ilia takes the form 

u — 1 + • • • 


G 0 Mx M ~ 1 


(4.15a, b) 


the resulting upper layer problem for (u, U) does not have a solution. The depth of the 
region IHb can now be inferred from equation (3.9). The dominant term on the right 
hand side of this equation is 0(x M ~ 1 ) so a balance with a comparable term on the left 
hand side can only be achieved if a + yi ~ X 2 ~ 3M . Thus 3/1 and a must be 0(x 2 ~ 3M ) and 
we therefore write 


2-3M 1 
Vi — J/ioz + 

„ _ „ J-3M , 

CL — CLqX + * 


(4.16a, b) 


The solution for (5, v) in IHb can then be developed in the form 


u — Nx 1 2M (yx 3M 2 ) + 


v = -{M -l)x M ~ 2 V- + 

It 


(4.17a, 6 ) 


Here N is an unknown constant to be determined. Finally the continuity of the normal 
velocity component at 3/1 requires that we expand b in the form 


h = bc,x 2 - 5M + • 


(4.17c) 
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The conditions (3.20), (3.21) for j — 1 then require that 



and 


\/ao + 2yio 

VGq 


Ny io, 


(M-l) 


J/io — ~bo + 


(3 M — 2)y/ap + 2yio M (ap + 2yio) 2 


2 a/ Gq 


+ 


VGo 


(4.18a, 6, c) 


It follows from (4.18a,b) that a 0 = — J/io so it remains for us to determine y 10 . The required 
equation follows from the zeroth order approximation to (3.9) which yields 


_ 9G 0 
Vl ° M 2 

and then N and b 0 follow from (4.17a,c) respectively Thus in the limit x —*■ oo with 
X ~ x M , M > 1/2 the locations t/i and y 2 have the asymptotic forms 


Vi = 


9G Q .2-3M , 
+ 


G 0 Mx m 

V2 = 2 f * ” 


(4.19a, b) 


and between these positions the mean velocity component in the x direction is given by 
the square root form (4.14). Thus for a curvature distribution \ = the initial and far 
downstream regions where a finite amplitude vortex exists are as indicated in Figure 3 for 
M > 1/2. The intermediate region can only be completed by a numerical procedure of the 
type we discuss in the next section. 


5. A numerical scheme for the determination of the free boundaries yi,y 2 


We shall now outline a scheme which we have used to solve the problem specified 
by (3.17) - (3.20). For convenience we drop the notation for u,v and assume that a 
solution of the problem is known for x < x. The scheme which we have used can be used 
to advance this solution downstream to x + e for sufficiently small values of e. In general 
the state for x < x must be calculated using the approach of §4 in a neighbourhood of 
where the vortices first become unstable. 


The first step in our calculation is to define a variable f by 


? = — J = 1 ,y < yi> 

Vj 


j = 2,y > y 2 , 
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so that (3.17) is now to be solved on (0, oo) in terms of f with boundary conditions at 0, 
and oo together with ‘jump’ conditions at f = 1. Thus (3.17) is now written in the form 

du _ 1 d 2 u v du y'j du 

U dx y'j d f 2 yy d$ + yy ^ d$ 


dv f du y'j du ) 

Vl l dx + f yy d$ j ’ 


(5.1<z, b) 


where j = 1 for f < 1 and j = 2 for f > 1. The required boundary conditions are 


u = v = 0, C = 0, w — » 1, f — >• oo 


(5.2) 


and the ‘discontinuous’ conditions at f = 1 may be written 


«± = 


\/o + 2yy 

%/^oX 


_ Vj 

V&oXy/a + 2yy ’ 


_ ayq + 2yy {q + 2yy>2x' 

± 2\/G r oX 6 %/GoX X 

Here the ± signs correspond to the limits f — ► 1 + , f — > 1- 
associated with the + sign and j = 1 with the negative sign, 
completed by the jump condition (3.9). 


(5.3a, 6, c) 

and the index j = 2 is 
Finally the system is 


We first advance the solution of (5.1a) for f in the range [0, 1]. This is done using the 
scheme 

i ■> {^n+1 2u n U n _i} = 


h 2 yj 


2 £Vn ( 1 . ^nj/y«n . . 

U " - 2^ { “" +1 _ + _ ^ < “" +1 _ 


(5.4) 


where h is the vertical grid spacing and a tilde denotes a quantity evaluated at the position 
x + e. The index n refers to a quantity at the grid point f = nh. In order to solve the 
tridiagonal system (5.4) we must make a guess for yy and set y'- = (yy — yy)/c. When 
solving (5.4) we satisfy the required condition on u at f = 0 and (5.3b) with j = 1. The 
continuity equation (5.1b) is then discretized as 


#n+i — £n-l 

2h 



CnVj U n +1 fin— 1 

Ti 2 h 


} 


(5.5) 


so that v can be determined at x + e for 0 < f < 1. However the equation for v is only 
of first order in £ so only the boundary condition at f = 0 can be satisfied during this 
procedure. Thus the solution of (5.1) for 0 < f < 1 has been calculated at x + e but, as 
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yet, (5.3a,c) with j = 1 have not been satisfied. However these conditions are now used to 
obtain an improved value of y x and a value for a by writing these conditions in the form 

2j/i = G 0 X«- ~ ~ a 


~ | l (a -a) \/d + Wi {d + 2 yi> 2 X / 

t ’- + 0 e eVG^x 


and iterating until a and y x converge. Here b is the current guess for b at x+ e. The scheme 
used to find u,v at x + e for 0 < f < 1 can be applied in a similar manner to the region 
f > 1. The u equation is solved subject to u — ► 1, f — * oo and (5.3b) with j = 2, whilst the 
v equation is solved subject to (5.2c) with j = 2. Finally (3.9), (5.3c) are written in the 
form 


K 


VZTWi 


\ /o + 2 j/2 1 

J 


x' 

12GqX 2 


{(a + 2y 2 ) 2 - (a + 2y x ) 2 } 


If A A ■, 

-| { — A 

\/G oX V® + 2y 2 y/a + 2y x 


2y2 = Gox«+ - o. 


The second of these equations determines a new value for y 2 and the first one then deter- 
mines a value for b. Thus we now have values for u, v, a, b, y x , y 2 at x + e. We then repeat 
the whole procedure using the new values of a, 6, y x ,y 2 obtained in the first iteration until 
converged values of these quantities are obtained. 


The above numerical scheme was found to converge for sufficiently small values of c 
the step length in the x direction. It was found that h = .1,1= .005 gave a stable scheme 
for the cases investigated and produced values for yi,y 2 and the other flow quantities 
correct to two decimal places. The first flow considered had Go = 5., x = v2x and the 
free boundary value problem of §3 was integrated for x > .5. At x = .5 the initial values 
of y x , y 2 etc. were calculated from the similarity solution of §3. In Figure 4 we have shown 
the values of y x , y 2 calculated using the scheme outlined in this section. In the same Figure 
we have indicated the values of yi,y 2 predicted by the similarity solution of §3. We see 
that there is excellent agreement between the results from the different solution methods. 

As an example of a non ‘self-similar’ Gortler vortex we considered the case x = 
2x,G 0 = 4.176. This curvature distribution has a basic state which is linearly unstable 
beyong x* = .5. The small (x — x*) asymptotic solution of §4 was used to generate initial 
values of yi,y 2 ,a,b and the velocity profiles at x = .51. The numerical scheme described 
in this section was then used to advance the solution beyond x = .51. The results obtained 
for y x and y 2 are shown in Figure 5 for .51 < x < 1.8. In this Figure we have also shown 
the corresponding results predicted by the small x — x* and large x asymptotic solutions of 
§4. We see that, apart from an apparently constant horizontal translation, the numerical 
scheme predicts values of y x and y 2 which rapidly approach their large x amplitude values. 
This translation is not unexpected since the large x solution of §4, to the order given, has 
an arbitrary origin for x. 
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In Figure 6 we have shown the streamwise velocity component produced by our scheme 
for the numerical solution of the free boundary problem in §3. The Blasius flow appropriate 
to this position and the large x asymptotic solution for u are also shown. It can be seen 
that at x = 1.8 the asymptotic and numerical solutions are virtually identical. Thus at 
larger values of x we can approximate u by the asymptotic solution of §4. In Figure 6d 
we have shown the result of making this assumption to calculate u at x = 7.5, 15. We see 
that by this step the boundary layer has been substantially thickened by the effect of the 
vortices. This is because the asymptotic solution has y 2 ~ X ~ x whilst the undisturbed 
boundary layer grows like x?. 

In Figures 7 and 8 we have shown the eigenfunctions Vq and Uq at the downstream 
locations x = .6, 1.0, 1.4, 1.8. The corresponding large x asymptotic solutions for Vq and 
Uq can be derived using the analysis of §4. For x = 1.8 such a calculation produces 
results virtually identical to the ones shown. In Figure 8 we see that when x increases the 
streamwise velocity component of theGortler vortex develops a shear layer near the lower 
boundary. The development of the shear layer is caused by the fact that the transition 
V 1 

layer lib approaches the wall when x — > oo. Thus Uq = increases rapidly from 

y/a + 2yi 

zero to an 0(1) value in the neighbourhood of y lt 

The implication of the above calculations is that we can reasonably expect that the 
large x solution of§4 will give accurate predictions for the vortex induced flow quantities at 
relatively small values of x. Thus in practical situations we might expect to obtain sensible 
results by using that approach rather than the numerical scheme for the free boundary 
problem. Moreover such an approach produces velocity profiles whose lower stability to 
Tollmien-Schlichting waves can be readily investigated. 

6. Conclusions 

It is perhaps useful at this stage to remind the reader that the small wavelength ap- 
proximation wehave used does not make our calculation physically unrealistic. This is 
because it is known experimentally that when Gortler vortices develop their wavelength 
is conserved as they move downstream. Thus for a growing boundary layer the effective 
wavenumber increases in the downstream direction so that a small wavelength approx- 
imation eventually becomes justifiable. We have no reason to suppose that in the non 
linear case the downstream position where the small wavelength results approach the 0(1) 
wavelength results will differ significantly from that for the linear case. 

The discussion in §3-5 has been concerned with flows for which x increases at least 
as quickly as ij when x increases. Otherwise the basic state will become linearly stable 
at a finite value of x and the structure we have found will terminate at some value of x. 
The termination of the finite interval of vortex activity is simply the ’mirror-image’ of the 
small (x — x*) solution of §4. This result can be confirmed from the weakly nonlinear 
theory of Hall (1982b) so that at some value x, say x + , r/j and y 2 will have the asymptotic 
form y + ± y ++ (x + — x)* where y + ,y ++ are constants. In Figure 9 we have shown the 
development of yi,y 2 for Go = 4.176 with x = 2x, x < 1, x = 4x 2 — 3x 3 , x > 1. We 
see that the region of vortex activity which begins at linear neutral position x = .5 stops 
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at x ~ 1.5. In the absence of a finite amplitude vortex beyond x = .5 this flow is linearly 
stable for x > 1.854. Thus the presence of the vortices in the range .5 < x < 1.5 produces 
a boundary layer which is stable in a regime where it would have been unstable if the 
vortices had not developed. 

However the effect of the vortices on the boundary layer does not end when the layer 
of vortex activity terminates. This is because the initial velocity distribution at x = x + 
for the subsequent boundary layer will in general be quite different from that appropriate 
to the undisturbed flow. Furthermore there is no reason to suppose that y + should equal 
y* the location of the vortices according to linear theory applied to the undisturbed state. 
Thus the decay of the vortices at x = x + does not in any sense allow the boundary layer 
to return to its undisturbed state. It follows that Gortler vortices might have a significant 
effect on separation subsequent to a region of concave curvature. Indeed it is known from 
the work of Hall and Bennett (1986) that triple-deck flows can support Taylor-Gortler 
vortices so the properties of these flows might also be significantly altered in the presence 
of vortices. 

In some flows it is possible that there will be several intervals in x where vortices can 
develop. We might expect that the steady boundary layer over a wavy wall might support 
vortices at regular intervals along the wall. In the intervals where vortices do not develop 
the basic state will in general be altered from its undisturbed state by the vortex activity 
in the previous undisturbed interval. 

The most surprizing feature of our calculation is that the fully nonlinear state driven 
by large amplitude Gortler vortices can be described in a relatively simple manner. The 
major effect of the vortices is to gradually expand into the boundary layer to give the 
mean flow there a simple square root profile. If the location increases faster than x? this 
layer thickens until it occupies the whole of the boundary layer apart from a thin layer at 
the wall. In addition the layer of vortex activity expands into the free stream and thus 
thickens the undisturbed boundary layer. 
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Figure 5. The development of y and y with x for the case x = 2x, g = 4.176 
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Figure 6d. The strearrwise velocity component corresponding to 
X = 2x, g Q = 4.176, x = 7.5, 15. 
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